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Abstract 

Properties of Braneworld models of dark energy are reviewed. Braneworld models 
admit the following interesting possibilities: (i) The effective equation of state can be 
w < —1 as well as w > — 1. In the former case the expansion of the universe is well 
behaved at all times and the universe does not run into a future 'Big Rip' singularity 
which is usually encountered by Phantom models, (ii) For a class of Braneworld 
models the acceleration of the universe can be a transient phenomenon. In this case 
the current acceleration of the universe is sandwiched between two matter dominated 
epochs. Such a braneworld does not have a horizon in contrast to LCDM and most 
Quintessence models, (iii) For a specific set of parameter values the universe can either 
originate from, or end its existence in a Quiescent singularity, at which the density, 
pressure and Hubble parameter remain finite, while the deceleration parameter and 
all invariants of the Riemann tensor diverge to infinity within a finite interval of cosmic 
time, (iv) Braneworld models of dark energy can loiter at high redshifts: 6 <, z < 40. 
The Hubble parameter decreases during the loitering epoch relative to its value in 
LCDM. As a result the age of the universe at loitering dramatically increases and 
this is expected to boost the formation of high redshift gravitationally bound systems 
such as 1O 9 M0 black holes at z ~ 6 and lower-mass black holes and/or Population III 
stars at z > 10, whose existence could be problematic within the LCDM scenario, (v) 
Braneworld models with a time-like extra dimension bounce at early times thereby 
avoiding the initial 'Big Bang singularity', (vi) Both Inflation and Dark Energy can 
be successfully unified within a single scheme (Quintessential Inflation). 

1 Introduction 

One of the most remarkable discoveries of the past decade is that the universe is accelerating. Evidence 
for an accelerating universe comes from observations of high redshift type la supernovae treated as 
standardized candles pQ and, more indirectly, by observations of the cosmic microwave background and 
galaxy clustering [21 E] . Perhaps the simplest explanation for acceleration is the presence of vacuum 
energy exhibiting itself as a small cosmological constant with equation of state p = —p = constant. 
However, its un-evolving nature implies that the cosmological constant must be set to an extremely small 
value in order to dominate the expansion dynamics of the universe at precisely the present epoch and 
this gives rise (according to one's perspective) either to an initial 'fine-tuning' problem or to a 'cosmic 
coincidence' problem. As a result several radically different alternative methods of generating 'dark 
energy' at a sufficiently late cosmological epoch have been suggested (see @J \5\ \G\ El E] for reviews on 
this subject). In this talk I will focus on one such approach which rests on the notion that space-time is 
higher-dimensional, and that our observable universe is a (3+l)-dimensional 'brane' which is embedded 
in a (4+l)-dimensional 'bulk' space-time. As we shall see, higher-dimensional braneworld models allow 
the expansion dynamics to be radically different from that predicted by conventional Einstein gravity in 
3+1 dimensions. Some cosmological 'surpises' which spring from Braneworld models include: 

• Both early and late time acceleration can be successfully unified within a single scheme (Quintessen- 
tial Inflation) in which the very same scalar field which drives Inflation at early times becomes 
Quintessence at late times. 

• The (effective) equation of state of dark energy in the braneworld scenario can be 'phantom-like' 
(to < —1) or 'Quintessence- like' (to > —1). (These two possibilities are essentially related to the 
two distinct ways in which the brane can be embedded in the bulk.) 
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• The acceleration of the universe can be a transient phenomenon: braneworld models accelerate 
during the present epoch but return to matter-dominated expansion at late times. 

• A class of braneworld models encounter a Quiescent future singularity, at which d — > constant, but 
a — > — oo. The surprising feature of this singularity is that while the Hubble parameter, density 
and pressure remain finite, the deceleration parameter and all curvature invariants diverge as the 
singularity is approached. 

• A spatially flat Braneworld can mimick a closed universe and loiter at large redshifts. 

• A braneworld embedded in a five-dimensional space in which the extra (bulk) dimension is time- 
like can bounce at early times, thereby generically avoiding the big bang singularity. Cyclic models 
of the universe with successive expansion-contraction cycles can be constructed based on such a 
bouncing braneworld. 

Let us now dwell a little on each of these cosmological properties ('surprises'). 



2 Quintessential Inflation on the Brane 

An intriguing issue in cosmology is that the universe appears to accelerate twice: once at the very 
beginning during Inflation and again about 10 billion years later, during the present epoch. Although 
most theoretical models assume that there is no real connection between the two epochs and that Inflation 
and Dark energy are distinct physical entities, it might well be that the two phenomena are in fact related, 
and that the same scalar field which initially drives inflation, later, when its density has been considerably 
reduced, plays the role of Quintessence. This notion of Quintessential Inflation was first explored in the 
context of the Einstein gravity in 3+1 dimensions by Peebles and Vilenkin The possibility that 

Braneworld models could provide a more efficient realisation of this scenario was suggested by Copeland, 
Liddlc and Lidsey ^2] and subsequently discussed in greater detail by several authors [T3l H^l IT51 ITfil ITT] . 

The main line of reasoning behind this approach is simple. In the Randall-Sundrum model |i] the 
modified Einstein equations on the brane contain high-energy corrections as well as the projection of the 
Weyl tensor from the bulk on to the brane. The Friedmann equation on the brane in this case becomes 

rr2 , « _ A b C (p+O-f , , 

H + ^-T + ^ + ^m^- (1) 

where k = 0, ±1, a is the brane tension, At, is the value of the cosmological constant in the five-dimensional 
'bulk' and C/a 4 is the 'dark radiation' term which describes the projected five-dimensional degrees of 
freedom onto the brane. It is easy to see from JIJ that 

# 2 (a-oo) = ^, (2) 

where 

Ars = T + 3^ (3) 

is the effective four dimensional cosmological constant in the Randall-Sundrum model. 

During Inflation the expansion of the universe is driven by a scalar field propagating on the brane 
with energy density and pressure given by 

p^Y + VW ' P* = Y~ v{(t>) ' (4) 

while the evolution equation for the scalar field is 

4>(t) + 3H(j>(t) + = . (5) 

Braneworld models add an interesting new dimension to the scalar-field dynamics. As demonstrated by 
P|l. for p> cr a quadratic density term appears in the modified Friedmann equation on the brane. This 
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radically alters the expansion dynamics at early epochs by speeding up the rate of expansion and the 
value of the Hubble parameter, which now becomes H oc p instead of the more conventional H cx J~p. 
Since the value of H affects the motion of tj> through the damping term l 3.ff</>' in (JSJ, the scalar field ex- 
periences significantly greater damping in braneworld cosmology than in conventional GR. Consequently, 
inflation on the brane can be realized by very steep potentials - precisely those that are used to describe 
Quintessence. The braneworld scenario therefore provides us with the opportunity to unify inflation and 
dark energy through the notion of quintessential inflation ^2 1121 1131 1141 1151 IT7| . An example of 
quintessential inflation is shown in Figure ^ 

An important property of Quintessential Inflation on the brane is that once the inflationary regime is 
over and braneworld corrections cease to play an important role, the extreme steepness of the potential 
causes the scalar field to plunge down its potential resulting in a 'kinetic regime' prior to reheating during 
which j»0 ~ p^. It is well known that the spectrum of relic gravity waves created quantum mechanicaly 
during Inflation is sensitive to and bears an imprint of the post-inflationary equation of state |18l I19| . 
The effect of the kinetic regime is to create a 'blue' gravity wave spectrum on short wavelength scales 
which is an important observational signature of Quintessential Inflation on the brane E| j see also 
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Figure 1: The post-inflationary density parameter Q is plotted for the scalar field (solid line), radiation (dashed line), and 
cold dark matter (dotted line) in the quintessential-inflationary model described by the potential 1221 V((j>) = Vb[cosh Acft— l] p 
with p = 0.2. Late-time oscillations of the scalar field ensure that the mean equation of state turns negative (w^) ~ —2/3, 
giving rise to the current epoch of cosmic acceleration with a(t) oc t 2 and present-day values Hrj^, ~ 0.7, f2rj m ~ 0.3. Note 
the existence of two epochs, early and late, when the scalar field dominates the matter density of the universe. Figure 
courtesy of 1171 . 



3 Braneworld models of Dark Energy with w < — 1 and w > — 1 

A radically different way of making the universe accelerate was suggested in [231 127| . The braneworld 
model which I now discuss presents a successful synthesis of the higher-dimensional ansatzes proposed 
by Randall and Sundrum [Sj and Dvali, Gabadadze, and Porrati [221, and is described by the action [2~4] 



S = AI S 



(11 - 2A b ) - 2 



bulk 
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L{h a b, 4>) ■ 
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Here, 1Z is the scalar curvature of the five-dimensional metric g a b in the bulk, and R is the scalar curvature 
of the induced metric h a b = g a b — n a nb on the brane, where n a is the vector field of the inner unit normal 
to the brane. The quantity K — K a i,h ab is the trace of the symmetric tensor of extrinsic curvature 
K a b = h c a \J c rib of the brane, and L(h a b, cf) denotes the Lagrangian density of the four-dimensional 
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matter fields 4> whose dynamics is restricted to the brane (the notation and conventions are those of |25p. 
Integrations over the bulk and brane are taken with the natural volume elements \J—g d 5 x and \J — h d 4 x, 
respectively. The constants M and m denote, respectively, the five-dimensional and four-dimensional 
Planck masses, Ab is the five-dimensional (bulk) cosmological constant, and a is the brane tension. 

An important difference between the action JBJ and that describing the Randall-Sundrum cosmology is 
the presence of the term m? J Rd i x in © . This term can be thought of as resulting from the backreaction 
of quantum fluctuations of matter fields residing on the brane, and this mechanism of generating the 
gravitational part of the action was proposed by Sakharov in a seminal paper in 1967 |2fi| . The effect of 
including such a term in the braneworld action was explored in |24l 1231 12 7| . 

Action (jSJ) leads to the following expression for the Hubble parameter on the brane for a spatially flat 
universe |27| : 




(7) 



where 



Note that the four-dimensional Planck mass m is related to the effective Newton's constant on the brane 
as m — 1/VSttG. The two Planck masses M and m define a new length scale t = 2m 2 /M 3 ~ cH 1 in a 
braneworld which begins to accelerate at the current epoch (23 E|- (Other applications of braneworlds 
to the late-time acceleration of the universe are discussed in |29l 1361 137)1 1311 ITTTTl 1341 I3"51 128j .) 

The two signs in JJJ) correspond to two branches of braneworld models and are related to the two 
different ways in which the brane can be embedded in the bulk. As shown in |27| . the '+' sign in Q 
corresponds to late time acceleration of the universe driven by dark energy with an 'effective' equation 
of state w > — 1 (we call this model BRANE2) whereas the '— ' sign is associated with phantom-like 
behaviour w < — 1 (this model is called BRANE1). 

Some limiting cases of (0) will be of interest to the reader: 

1. For m = equation (JJJ) reduces to a Randall-Sundrum universe ^TJ described earlier in 0J. 

2. In the other extreme case when M = 0, extra-dimensional effects become unimportant, and J7J 
reduces to the LCDM model 

H 2 (a) = ^+B. (9) 

ar 

3. Finally, when Ab = and a = 0, Q reduces to the DGP braneworld |23|. 

Consider now the Braneworld with the ' — ' sign on the right-hand side of Q, i.e. BRANE 1. In this 
case the Hubble parameter on the brane can be rewritten as: 

H 2 (a) = — + A eff . 
or 

The term A e g is composed of two terms, namely, a constant A-term and a 'screening term': 



A eff H£ + f)-f f + ^ + (10) 



4 a- (ii) 

A Screening term 

Since the screening term decreases with time, the value of the effective cosmological constant 
A c ff increases. Therefore the Braneworld behaves just like a Phantom model (w < —1) but without 
Phantom's problems: no violation of the weak energy condition and no future singularity. Indeed, from 
ifTUjl its clear that the universe evolves to ACDM in the future fsee also 36. 40. 41. 371 1361 13131 l4Hl4Hl4l| 
113 ESI 113 SHI SHI EOI EH E3 E31 for discussions of related issues) . 
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The fact that the Braneworld model J7J) can give rise to Phantom-like behaviour can also be seen if 
we rewrite Q in terms of the cosmological redshift 'z' so that |77) 



H 2 (z) 



Hi 



2^h V«m(l+^) 3 + n a + n e + Q Ab , (12) 
where we have set the dark radiation term in (0 to zero (C = 0) and 

o _ 1 o _ P° m o _ Z. o _ Ab f-\o\ 
1 ~ PHI ' m " 3m*ff 2 ' ~ 3m 2 ^ 2 ' Ab ~ 6H$ ' 1 ' 

The quantities f2j satisfy the constraint equation 



where the ± sign in l|14(l is commensurate with that in ljl2|l . For future reference we also note the 
cosmological density associated with the four-dimensional cosmological constant (J2J) in Randall-Sundram 
cosmology 

n—^-^-^- (15 » 

The expression for the current value of the effective equation of state is 



w 2g - 1 = n m / 

w ° 3(1 -a„) l-fim V n m + n,, + n f + n Ab ' 1 j 

and we immediately find that wq < — 1 when we take the lower sign in (|lfj|) . which corresponds to choosing 
one of two possible embeddings of this braneworld in the higher dimensional bulk (this is BRANE1). (The 
second choice of embedding (BRANE2) gives wq > — 1.) It is important to note that all BRANE1 models 
have wo < —1 and w(z) ~ —0.5 at z 1. They therefore successfully cross the 'great divide' at w = — 1 
(see page 26 of j2ZJ for a discussion of this issue). 



4 Transient Dark Energy 

Braneworld cosmology has the interesting property that it allows the acceleration of the universe to be 
a transient phenomenon. (In these models were referred to as 'Disappearing Dark Energy' (DDE).) 
Indeed, by setting the effective four dimensional cosmological constant (j2J to zero so that fJpts = in 
HHJ we get 

n a = ±2^uiZ. (17) 

Selecting the lower sign in <|17(l and substituting Q a = — 2-^/^12 A b in i|12|) we find that 2 



H(z 



Hn 



2^n e \M-v^± vtot-vn 



(18) 



which vanishes for the upper sign when < ^A b - This implies that, since the effective four dimensional 
cosmological constant vanishes, the universe reverts to matter dominated expansion in the future: 
H 2 (z — > —1) ~ fi m (l + z) 3 — > 0. The current acceleration of the universe is therefore a transient phase 
in this model, as illustrated in figure 

It is well known that an eternally accelerating universe is endowed with a cosmological event horizon 
which prevents the construction of a conventional S-matrix describing particle interactions within the 
framework of string or M-theory |54l . A transiently accelerating braneworld may therefore help reconcile 
string/M-theory with accelerating cosmology, which is an attractive feature of this braneworld cosmology. 



2 Note that the case = 2^yfl(Q Ab leads to transient acceleration with Q, m > 1 which is unrealistic 1271 . 
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Figure 2: The deceleration parameter g = —d/aH 2 is plotted as a function of redshift for the BRANE2 model 
with the Randall-Sundrum constraint = — 2^jQ,iVl^ h . The vertical (dot-dashed) line at z = marks the 
present epoch, while the horizontal (dashed) line at q — corresponds to a Milne universe [a(t) oc t] which neither 
accelerates nor decelerates. The case q = 0.5 corresponds to matter dominated expansion. Note that the universe 
ceases to accelerate and becomes matter dominated in the past as well as in the future. Figure courtesy of |27|. 

5 Quiescent Future Singularity 

As remarked earlier, an interesting property of the Braneworld J7J) is that the late-time expansion of the 
universe can culminate in a 'Quiescent' singularity at which the density, pressure and Hubble parameter 
remain finite while the deceleration parameter and geometrical invariants constructed from the Riemann 
tensor diverge [55] . 

In order to appreciate the origin of such unusual future singularities in the braneworld, consider again 
the expansion law l|12f) . For simplicity, we shall only discuss the solution corresponding to the '+' sign 
in O (called BRANE2 in 

A necessary and sufficient condition for the existence of a Quiescent singularity is for the inequality 



to be satisfied [55]. In this case the expression under the square root of I112[I becomes zero at a suitably 
late time, and the cosmological solution cannot be extended beyond this time. It can be shown that the 
scale factor a(t) and its first time derivative remain finite, while all the higher time derivatives of a(t) 
tend to infinity as the singularity is approached. 

The limiting redshift, z s = ao/a(z s ) — 1, at which the braneworld becomes singular is given by 



The time of occurance of the singularity (measured from the present moment) can easily be determined 
from 



n a + n e + n Ab < o , 



(19) 
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Figure 3: The deceleration parameter (solid line) is shown for a braneworld model with 57 m = 0.2, Qi = 0.4, 
Q.A = fi K = 0, and Q. a determined from 114H . We find that q(z) — ► 0.5 for z 2> 1 while g(z) — > 00 as 2 — > 
—0.312779... Currently go < 0, which indicates that the universe is accelerating. Also shown is the dimensionless 
Hubble parameter h(z) = 0.1 x H(z)/Ho (dashed line) which remains finite in this model. The vertical line at 
z — shows the present epoch. Figure courtesy of |55|. 



where H(z) is given by l|12|) . In Fig. [31 we show a specific braneworld model having f2 m = 0.2, fig = 0.4, 
r^Ab = = 0. In keeping with observations of high redshift supernovae our model universe is currently 
accelerating pQ, but will become singular at z s ~ —0.3 a(z s ) ~ 1.4 x ao, i-e. after T s ~ 4.5 /i^ q Gyr 
(frioo = i?o/100 km/sec/Mpc). Figure |5] demonstrates that the deceleration parameter becomes singular 
as z s is approached: 

a iJ' 

g = - T77I = — (l + z)-l; lim g(z) -> 00 , (22) 



aiJ 2 i? 

while the Hubble parameter remains finite: 



H 2 {z s ) 
Hi 



n e - r» Ab . (23) 



The possible presence of such 'Quiescent' singularity in Braneworld theory can also be seen from the 
following consideration. As shown in [37j, the action © leads to the bulk being described by the usual 
Einstein equation with cosmological constant: 

Gab + A b5ah = , (24) 

while the field equation on the brane is 

m 2 G ab + ah ab = r ab + M 3 (K ab - h ab K) . (25) 

Here, T ab is the stress-energy tensor on the brane. A divergent form for the extrinsic curvature (underlined 
term) caused by a singular embedding of the brane in the bulk will, through l|25|) , lead to a singular value 
for the Einstein tensor G ab even though all components of the stress-energy tensor on the brane remain 
finite ! 

More generally it is well known that Quiescent singularites occur when the original equations of motion 
are non-linear with respect to the highest derivative. They have earlier been discussed in the context of 
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Einstein gravity with the conformal anomaly 56 . (This result is not surprising in view of the formal 
similarity between braneworld theory and GR-based models with the conformal anomaly, discussed in 
|27|.) 'Determinant singularities' having a similar structure and properties are known to arise in the 
anisotropic Bianchi I model containing a dilaton coupled to a Gauss-Bonnet term in the action |57) and 
also in other cosmological models [SSJ . For instance [SS] refer to singularities in which the deceleration 
parameter tends to infinity as the 'Big Brake' while in |5H| they are called 'sudden' singularities (see also 

EDI). 

We should also emphasise that both the past and future quiescent singularities occur for a wide range 
in parameter space and might provide an interesting alternative to the 'big bang'/'big crunch' singularities 
of general relativity. 

6 Loitering Braneworld Models 




Figure 4: Dynamical possibilities in a closed FRW universe are essentially the same as in the spatially flat brane 
and are summarised above. Figure courtesy of pT|. 

An interesting aspect of the Braneworld models J7J) is that they allow for the existence of an early 
loitering epoch |61| . Loitering is characterized by the fact that the Hubble parameter dips in value over 
a narrow redshift range referred to as the 'loitering epoch'. During loitering density perturbations are 
expected to grow rapidly. In addition, since the expansion of the universe slows down, its age near 
loitering dramatically increases. An early epoch of loitering is expected to boost the formation of high 
redshift gravitationally bound systems such as 10 9 M Q black holes at z ~ 6 and lower-mass black holes 
and/or Population III stars at z > 10, whose existence could be problematic within the LCDM scenario. 

To demonstrate the possibility of loitering note that for large values of the 'dark radiation' term C/a 4 , 
when C < and £ 2 \C\/a 4 ^> 1, equation Q acquires the form 

H\a)*^ + B± 2 -^. (26) 



Equation l|26|l bears a formal resemblance to the Hubble parameter in standard GR 



■> ... 3 n z 



Remarkably, the role of the spatial curvature in l|27|l is played by the dark-radiation term in i|26|) : con- 
sequently, a spatially open universe is mimicked by the BRANE2 model whereas a closed universe is 
mimicked by BRANE1. 




Figure 5: The Hubble parameter for three different universes all of which loiter at zi - lt ~ 18. Parameter values 
are f2 m = 0.3, Q.c = 8.0, fit = 3.0, and fiA b /10 5 = 6, 4.5, 3.4 (solid lines, from top to bottom). The left panel 
shows the Hubble parameter with respect to the LCDM value while, in the right panel, the LCDM (dashed) and 
loitering (solid) Hubble parameters are shown separately. Figure courtesy of |61| . 



Loitering solutions in standard cosmology were first found for a closed FRW model with a cosmological 
constant in | 62| . Loitering in more general situations was discussed in |63l 1641 1651 161j . An intriguing 
aspect of braneworld models is that, since a spatially flat braneworld can behave dynamically just like 
a closed universe, loitering may be realised within a spatially flat setting and therefore be consistent 
both with Inflationary cosmology and recent CMB results 2 . Indeed, loitering solutions to 10 & (|26[) 
can easily be found by requiring that the loitering condition dH/dz = is satisfied. This allows us to 
determine the 'loitering redshift' jBJ 

4 yn^r e 

1+zioit "3^r- (28) 

From this expression we find that the universe will loiter at a large redshift z\ \t 3> 1 provided flc^i 3> ^m- 
Since O^, <C 1, loitering at large redshifts is not difficult to achieve. 

One should note here that a closed FRW universe with a cosmological constant can loiter only at 
rather small redshifts: z\ i t < 2 for Q ro > 0.1 [S5| . Therefore Braneworld cosmology is endowed with two 
fundamentally new attributes: (i) it permits loitering in a spatially flat universe and (ii) it allows the 
loitering redshift to be large. Neither of these possibilities is allowed in standard Einsteinian cosmology. 

The Hubble parameter at loitering is given by the approximate expression |oT| 

#W)„ n _ 32(^)3/2 

H* 27 Ql • [ ' 
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Examples of a loitering model are shown in Fig. [S] where the Hubble parameter of a universe which 
loitered at z ~ 18 is plotted against the redshift. The right hand panel of figure [S] illustrates the fact that 
the loitering universe can show a variety of interesting behaviour: (i) top curve, H (z) is monotonically 
increasing and H'(z) ~ constant in the loitering interval; (ii) middle curve, H(z) appears to have an 
inflexion point (H' ~ 0, H" ~ 0) during loitering; (iii) lower curve, H{z) has both a maximum and a 
mininimum, the latter occuring in the loitering regime. 
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Figure 6: The age of three loitering models is shown relative to the age in LCDM (model parameters are the 
same as in Fig. O . Note that the age of the universe near loitering (zioit ~ 18) is significantly greater than that in 
LCDM although, at the present epoch, the difference in ages between the two models is relatively small. Figure 
courtesy of 

A loitering universe could have several important cosmological consequences, one of them being that 
the age of the universe during loitering increases, as shown in Fig. |SJ The reason for the age increase 
rests in the expression 

*> = f (TTW < 30 > 

which shows that a lower value of H{z) close to loitering will boost the age of the universe at that epoch. 
An important consequence of having a larger age of the universe at z ~ 20 (or so) is that astrophys- 
ical processes at these redshifts have more time in which to develop. This is especially important for 
gravitational instability which forms gravitationally bound systems from the extremely tiny fluctuations 
existing at the epoch of last scattering. Thus, an early loitering epoch may be conducive to the formation 
of Population III stars and low-mass black holes at z ~ 17 and also of ~ 10 9 M Q black holes at lower 
redshifts (z ~ 6) whose existence could be problematic within the LCDM scenario [HSj- (Note that age 
of a LCDM universe at z > 1 is t(z) ~ (2/3iJ \/TQ(l + z)~ 3 ^ = 5.38 x 10 8 (1 + z/10)- 3 / 2 years for 
il m = 0.3, h = 0.7.) 

Another important property of a loitering universe is that it can alter the reionization properties of 
the inter-galactic medium at moderate redshifts z < bolter. One should note in this context that a major 
surprise emerging from the WMAP experiment was the discovery of a large optical depth to reionization 
t = 0.17 ± 0.06 0, which, when translated within the framework of LCDM (assuming instantaneous 
reionization), implies a rather early epoch for reionization z re ion — 17 ± 5. 

In order to appreciate how loitering might alter these conclusions consider the following expression 
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for the optical depth to a redshift z re i on 

T(z relon ) = c / — — — — (31) 
Jo ( l + z ) H ( z ) 

where n e is the electron density and o~t is the Thompson cross-section describing scattering between elec- 
trons and CMB photons. During loitering H(z) drops below its value in LCDM, therefore z re ion| loitering < 
Zreion|LCDM- For instance the redshift of reionization drops to z re i on < 12 for the loitering models shown 
in figure [SJ Loitering decreases the redshift of reionization and increases the age of the universe thereby 
alleviating the existing tension between the high redshift universe and dark energy cosmology |61) . 

6.1 The effective equation of state in a loitering Braneworld. 

The deceleration parameter q and the effective equation of state w in our loitering model are given by 
the expressions 

3 [1 - U m (z)\ 

where H(z) is determined from (|12fl . The current values of these quantities are 



qo = 2 S2r 



- 1 



where Y,i &i = + + + ^<x + ^A b = (V^h + V 1 + + ^c) • From Eq. we find that 
wo < —1 if f^c > 0; in other words, our loitering universe has a phantom- like effective equation of state 
(EOS). (In particular, for the loitering models shown in Fig. we have u>o = —1.035, —1.04, —1.047 
(top to bottom), all of which are in excellent agreement with recent observations |^7||68 .) 

An interesting consequence of the loitering braneworld is that the matter density fl m = 8irGp m /3H 2 
exceeds unity at some time in the past. This follows immediately from the fact that, since the value of 
H(z) in the loitering braneworld model is smaller than its counterpart in LCDM, the value of Q m (z) is 
larger than its counterpart in LCDM. One important consequence of this behaviour is that, as expected 
from (J2SJ, the effective equation of state (EOS) blows up precisely when fi m = 1. 

In Fig. [7] we show that, in contrast to the singular behaviour of the EOS, the deceleration parameter 
remains finite and well behaved even as w — > oo. Note that the finite behaviour of q(z) when Q nl (z) = 1 
reflects the fact that the EOS for the braneworld is an effective quantity and not a real physical property 
of the theory. This is an important point with broader ramifications since it concerns dark energy models 
constructed from modifications to the 'gravity sector' of the theory of which Braneworld models are but 
one example (others being: scalar-tensor gravity, Cardassian approximation, R+ f(R) gravity, k-essence, 
to name but a few). Let me therefore dwell on this point in some more detail. 

Within the framework of general relativity the rate of expansion of a FRW universe and its acceleration 
are described by the pair of equations 

o 8irG sr-^ k 

i 

d 4ttG 



a 



+ 3ft), (34) 



where the summation is over all matter fields contributing to the dynamics of the universe. Within this 
framework, and in the absence of spatial curvature, the energy density and pressure of dark energy can 
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Figure 7: The effective equation of state of dark energy (solid) and the deceleration parameter (dashed) are 
shown for a universe which loitered at z ~ 18. Note that the effective equation of state of dark energy becomes 
infinite at low redshifts when fi m (z) = 1. However, this behaviour is not reflected in the deceleration parameter, 
which becomes large only near the loitering redshift when H ~ 0. Figure courtesy of |ril| . 



be denned as: 

3H 2 

PDE — Political — Pm = „ (1 — ^m) , 

07TG 

H 2 ( l\ 

Pde = ^(ff-gj > (35) 

where Pcritical = 3H 2 /8irG is the critical density associated with a FRW universe and q = —a/aH 2 is the 
deceleration parameter. An important consequence of using (|34|1 & (|35|1 is that the ratio wde = Pde/ Pde 
can be expressed in terms of the deceleration parameter 

, , 2q(x) - 1 (2x/3) d InH I dx - 1 
wcs(x)= qK ' = \ 11 ' 3- , x = l + z. 36 

3(i-o m (x)) i - (H Q /Hyn m x 6 

i.e. we recover equation (|32f) which had earlier been used to determine the equation of state in a loitering 
universe. Since this derivation of w e e(z) uses the general relativistic formulation it can be expected to 
recover standard results for GR-based dark energy models such as LCDM and quintessence. However the 
same cannot be said of braneworld models since the Hubble parameter for the latter contains interaction 
terms between matter and dark energy (see for instance equation (|12|l ) and therefore does not subscribe 
to the general relativistic format 13411 . "One can however extend the above definition of w e g to non- 
Einstcinian theories by defining dark energy density to be the remainder term after one subtracts the 
matter density from the critical density in the Einstein equations. It should be emphasised that, according 
to this prescription all interaction terms between matter and dark energy (such terms arise in scalar-tensor 
and brane models) are attributed solely to dark energy. Therefore w e s(z) defined according to (|36|l is an 
effective equation of state in these models and not a fundamental physical entity (as it is in LCDM, for 
instance)" [021 ■ The fact that the EOS is indeed an effective equation of state and not a fundamental 
physical quantity is illustrated by figure which shows that w e s(z) diverges when Vt m {z) ~ 1 and that 
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this divergence is not reflected in other characteristics of expansion such as the deceleration parameter. 

3 



6.2 The Statefinder diagnostic and Dark Energy 

As demonstrated above, the equation of state is not a fundamental physical quantity for dark energy 
models based either on non-Einsteinian gravity (scalar-tensor theories, Cardassian approximation, k- 
essence, R + f(R) gravity, etc.) or on Einstein gravity in more than 3+1 dimensions, or indeed in 
cosmological models with explicit interaction terms between dark matter and dark energy. In all these 
cases it is useful to supplement the effective equation of state by other geometrical variables, which are 
'model independent', in the sense that they do not depend upon an underlying theory of gravity in an 
explicit way, but can be constructed solely from the expansion factor and its derivatives, and therefore 
would be expected to apply to all metric theories of gravity defined on a space-time manifold. In all such 
theories the universe can be expected to be characterised by the following general form for the expansion 
factor: 

o(t) = a(t ) + a\ Q (t - t ) + °^-{t - t ) 2 + ^i(f - t ) 3 + ... . (37) 

Usually dark energy models such as quiessence, quintessence, k-essence, braneworld models, Chaplygin 
gas etc. give rise to families of curves a(t) having quite different properties. Since we know that the 
acceleration of the universe is a fairly recent phenomenon we can, in principle, confine our attention to 
small values of |i — to\ in l|3"?Jl. 

Accordingly, as dicussed in [7DHI3], a new diagnostic of dark energy called 'statefinder' can be con- 
structed using both the second and third derivatives of the expansion factor (see also |721 1731 17^] ). The 
statefinder pair {r, s}, defines two new cosmological parameters (in addition to H and q): 

o - 9iu^ ,„ . 3^ w ,„„. 

= i + — n x (i + w)--n x - , (38) 



aH 3 2 v ' 2 A H 

r — 1 1 w , , 

3 55 %^i72T = 1 + --3^- (39) 

In figure[S]we show the statefinder pair for several dark energy models including LCDM, quintessence 
and the Chaplygin gas. For the sake of completeness we say a few words about each of these models 
below. (Note that the ensuing discussion is largely based on |7(J1 169j and the reader should refer to these 
papers for more details.) 

• In Quiessence the equation of state of dark energy is a constant so that 

H(z) = H [n m (l + zf + fl x {l + z) 3 ( 1+w A 1/2 . (40) 

where w — constant. Important examples of quiessence include: LCDM (w — —1), a network 
of non-interacting cosmic strings (w = —1/3), domain walls (w = —2/3) and phantom models 
(w < —1). Quiessence in a FRW universe can also be produced by a scalar field (quintessence) 

-2(l + m) 

which has the potential V((j)) oc sinh «» (C(j) + D), with appropriately chosen values of C and D 

• The density and pressure of Quintessence are given by 

P^l^ + V^), p^ = ^ 2 -V(0) , (41) 
while evolution of the Quintessence field <f> is governed by the equation of motion 

4> + 3H<j) + %■ = , (42) 



3 It should also be stressed that the propagation velocity of small inhomogeneities in dark energy is generically neither 
\/ w DE) nor \J <h>D E I dpDE ■ Therefore although w(z) is an important physical quantity it does not provide us with an 
exhaustive description of dark energy and its use as a diagnostic should be treated with some caution, as emphasised in 
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Figure 8: The left panel (a) shows the time evolution of the statefinder pair {r,s} for quintessence models and 
the Chaplygin gas. Quintessence models lie to the right of the LCDM fixed point (r = 1, s = 0) (solid lines 
represent tracker potentials V = Vo/(f> a , dot-dashed lines representing quiessence with constant equation of state 
to). For quiessence models, s remains constant at 1 + w while r declines asymptotically to 1 + |iu(l + w). For 
tracker models, s monotonically decreases to zero whereas r first decreases from unity to a minimum value, then 
rises to unity. These models tend to approach the LCDM fixed point (r = l,s = 0) from the right at t — > oo. 
Chaplygin gas models (solid lines) lie to the left of the LCDM fixed point. The right panel (b) shows the time 
evolution of the pair {r, q}, where q is the deceleration parameter. The solid line, which corresponds to the time 
evolution of the LCDM model, divides the r — q plane into two halves. The upper half is occupied by Chaplygin 
gas models, while the lower half contains quintessence models. All models diverge at the same point in the past 
(r = 1, q= 0.5) which corresponds to a matter dominated universe (SCDM), and converge to the same point in 
the future (r = 1, q = —1) which corresponds to the steady state model (SS) - the de Sitter expansion (LCDM 
— > SS as t — > oo and Q, m — » 0). The dark dots on the curves show current values {ro,so} (left) and {ra,qo} 
(right) for different dark energy models. In all models, Q m = 0.3 at the current epoch. In both panels quiessence 
is shown as dot-dashed while dashed lines mark envelopes for Chaplygin gas (upper) and quintessence (lower). 
Figure courtesy of |69|. 



where 



8ttG 



p 0m (l + z) 3 + -4> 2 + V(q>) 



(43) 



It is clear from (I41|l that w < —1/3 provided cb 2 < V(<p). Models with this property can lead 
to an accelerating universe at late times. An important subclass of Quintessence models displays 
the so-called 'tracker' behaviour during which the ratio of the scalar field energy density to that 
of the matter/radiation background changes very slowly over a substantial period of time. Models 
belonging to this class satisfy V"V / (V) 2 > 1 and approach a common evolutionary 'tracker path' 
from a wide range of initial conditions. 

• An interesting alternative form of dark energy is provided by the Chaplygin gas [731 [ZH] which 
obeys the equation of state 

Vc = -A/p c . (44) 
The energy density of the Chaplygin gas evolves according to 



Pc = ^A + B{l + zf 



(45) 



14 




Figure 9: Trajectories in the statefinder plane {r, q} for the Braneworld models discussed in 1121 . BRANE1 
models have w < — 1 generically, whereas BRANE2 models have w > — 1. The closed loop represents 'Disappearing 
Dark Energy' (DDE) (discussed in section|IJ for which the acceleration of the universe is a transient phenomenon. 
The thick solid curve in BRANE2 corresponds to the DGP model |23| . Dark dots indicate the current value of 
{r, q} for the models. All models are in reasonable agreement with current supernova data. One can see that the 
statefinder plane {r, q} easily helps distinguish the three Braneworld models from each other. Figure courtesy of 

m- 

from where we see that p c — > \J~A as z — > — 1 and p c — > y/B(l + z) 3 as z > 1. Thus, the Chaplygin 
gas behaves like pressureless dust at early times and like a cosmological constant during very late 
times. 

The Hubble parameter for a universe containing cold dark matter and the Chaplygin gas is given 

by 

I "I 1/2 

H(z)=H a rUl + z) 3 + ^^ + (l + z)6 . ( 4(>) 
where k = po m /yB. It is easy to see from (|46|) that 

POm i \ 

K = (Z — > OOj . 

Pc 



(47) 



Thus, k defines the ratio between CDM and the Chaplygin gas energy densities at the commence- 
ment of the matter-dominated stage. 

For all of the above models the statefinders r and s can be easily expressed in terms of the Hubble 
parameter H(z) and its derivatives as follows: 



H' 

r[x) = 1 — 2— x 




s(x) 



r(x) — 1 
3(q(x) - 1/2) ' 



(48) 
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where x = 1 + z and H is given by ((4011 . 143|) . (|46|) . 1(12(1 for the different dark energy models. The results, 
shown in figures|S]&|ni demonstrate that the Statefinder diagnostic can successfully differentiate between 
dark energy models as diverse as LCDM, Quintessence, the Chaplgyn Gas and Braneworlds (see [77] for 
other recent applications of the Statefinder diagnostic). 



7 Bouncing Braneworlds 

A remarkable feature of Braneworld cosmology is that the expansion dynamics of the early universe can 
be non-singular. Indeed, as shown in |78| bouncing Braneworld models can be constructed, based upon 
a Randall-Sundrum type action, but in which the extra dimension is time-like. In this case one starts 
from the action [7S| 



f (11 - 2A) - 2e f K - f 2a + f L(h a 

J bulk J branc J J brane J brane 



(49) 



where the parameter e = 1 if the signature of the bulk space is Lorentzian, so that the extra dimension 
is spacelike, and e = —1 if the signature is (— , — , +, +, +), so that the extra dimension is timelike. The 
evolution equations resulting from (|49|l have the form |78| 

2 _ Aeff 8irG NP ep 2 C 

Where Gn = ecr/127rAf 6 is the effective gravitational constant, and A e g = A/2 + eer 2 /3M 6 is the effective 
cosmological constant on the brane. For e = 1, Eq. i|50|) reduces to the well-known Randall-Sundrum 
form (we ignore n/a 2 ,A e g and C/et 4 which are unimportant for our purposes) 

On the other hand if e = — 1 then we obtain 

2 87tGn ( , p 



where the brane tension a < 0. (The two braneworld models 1)51(1 and l|52() are dual, as demonstrated in 
|79j.) Equation 152(1 clearly demonstrates that H ~ when pbouncc = 2 j crj i.e. the universe bounces when 
the energy density of matter becomes sufficiently large. It is important to note that the singularity-free 
nature of the early universe is a generic outcome of this theory and does not depend upon whether or 
not matter violates the energy conditions. An example of a bouncing braneworld universe is shown in 
figure 1(10(1 below. 

An early time bounce, such as 1(52(1 . can be used to construct cyclic models of the universe (8011811 1^1 
I88( . Consider for instance the mechanism proposed in |80j in which Phantom dark energy is presumed 
to exist in addition to matter and radiation so that 

P = Pr + Prn + PP (53) 

where p r cx a~ 4 , p m oc a~ 3 , pp cx a~ 3( - 1+,i '-' and w < —1 in the case of Phantom leading to 

p=4 + 4 + C« 3(|1+H) ■ (54) 

It is clear from 1(54(1 that p grows at small as well as large values of a(t). Indeed within one cycle the 
value of H will pass through zero twice: (i) at early times when the bounce in (I52() is brought about due 
to large values of the radiation density and (ii) during late times, when the large value of the Phantom 
energy leads to H = in 1(52(1 and initiates the universe's recollapse. 

Similarly, it is easy to show that a bouncing spatially closed universe with matter satisfying p + 3p > 
will also be 'cyclic' in the sense that it will pass through an infinite number of nonsingular expanding- 
collapsing epochs (811 152) . As demonstrated in |5J , a massive scalar field in such a universe usually 
leads to an increase in the amplitude of consecutive expansion cycles and to a gradual amelioration of 
the flatness problem; see figure ITTI 
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Figure 10: A bouncing radiation-dominated braneworld (rj = J dt/a is the conformal time). Figure courtesy of 

m 



8 Conclusions 

Braneworld models hold interesting consequences both for the early as well as the late-time evolution of 
the universe. At early times, the increased rate of expansion in RSII type braneworlds allows scalar fields 
with steep potentials to play the dual role of being the inflaton as well as Quintessence. 

Braneworld models in which the RSII action is supplemented by an induced gravity term on the brane 
have important consequences for the late-time evolution of the universe. For instance: (i) The effective 
equation of state of dark energy in this scenario can be w < — 1 as well as w > — 1. (ii) The current 
acceleration of the universe can be a transient phenomenon, (iii) The braneworld universe can end its 
existence in a Quiescent singularity at which the density, pressure and Hubble parameter remain finite, 
while the deceleration parameter and all invariants of the Riemann tensor diverge, (iv) Braneworld models 
of dark energy can loiter at high redshifts: 6 <, z < 40. The Hubble parameter decreases during the 
loitering epoch relative to its value in LCDM. As a result the age of the universe at loitering dramatically 
increases and this is expected to boost the formation of high redshift gravitationally bound systems such 
as 10 9 M e black holes at z ~ 6 and lower-mass black holes and/or Population III stars at z > 10, whose 
existence could be problematic within the LCDM scenario, (v) In addition to the above, an RSII type 
Braneworld in which the extra dimension is time-like (instead of space-like) avoids the initial Big Bang 
singularity by bouncing at early times. This property (which holds even if matter satisfies all the energy 
conditions) has been used to construct cyclic models of a singularity free universe. 
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